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Abstract 

We derive a novel and rigorous correction to the classical Reynolds lubrication approximation for fluids with 
viscosity depending upon the pressure. Our analysis shows that the pressure dependence of viscosity leads 
to additional nonlinear terms related to the shear-rate and arising from a non negligible cross-film pressure. 
We present a numerical comparison between the classical Reynolds equation and our modified equation and 
conclude that the modified equation leads to the prediction of higher pressures and viscosities in the flow 
domain. 
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1. Introduction 

The Reynolds equation [1], which is an approximation of the classical Navier-Stokes equations, describes 
reasonably well the flow of a large class of fluids whose viscosities can be assumed to be independent of the 
pressure in many lubrication problems. However, there is a clear and incontestable evidence, that is well- 
documented, that attests to the fact that the viscosity varies with pressure, especially so in several problems 
concerning thin film lubrication, cf. [2, 3, 4]. Experiments have shown that in high pressure regimes pressure 
variations can significantly change the viscosity of certain lubricants and in areas such as elastohydrodynamic 
lubrication the classical isoviscous lubrication theory is noticeably incapable of explaining the existence of 
continuous lubricant films, for example, in rolling-contact bearings, cf. [4J. 

The traditional correction to the Reynolds approximation in the piezoviscous regime is based on a 
rather heuristic assumption that it suffices to replace the constant viscosity in the Reynolds equation by 
a suitable viscosity-pressure relationship, cf. [5]. This approach becomes questionable, however, in high 
pressure regimes because of the possible change of type (loss of ellipticity) of the equations and the potential 
existence of cross-film pressure gradient, see the discussion in [6]. In fact, for a Reynolds type approximation 
to be valid in this regime it ought be derived from the full balance of linear momentum equations governing 
the flow of incompressible fluids with pressure-dependent viscosities. Now, the mathematical theory for the 
equations that govern the flows of fluids with a pressure-dependent viscosity has been advancing in leaps 
and bounds over the last few decades, see [7, 8, 9, 10, 11, 12, 13, 14, 15], but the only rigorous attempt to 
derive a modified Reynolds equation for elastohydrodynamic lubrication based upon these equations seems 
to have been carried out by Rajagopal and Szeri [6], see also [16] for further applications of their model. 

In this paper we propose a new modified Reynolds equation for hydrodynamic lubrication in high pressure 
regimes. Our approach is built upon an asymptotic expansion of the non-dimensional velocity and pressure 
fields in terms of a small dimensionless parameter e, related to the film thickness, and on a systematic 
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analysis of the simplified set of equations obtained at different orders of e, see [17] for a similar approach in 
the isoviscous case. Assuming that the dimensionless pressure-viscosity coefficient is of order e, we show that 
the pressure distribution is governed by a Reynolds equation modified by a term depending, in particular, 
on the square of the shear rate. In [6], the authors assumed, for simplicity, that the cross-film pressure 
vanishes and derived a modified Reynolds equation, similar to ours, but depending on the elongation rate. 
Our computations show that it is exactly the (lower-order) cross-film pressure which is responsible for the 
new modified term in the Reynolds equation. We also prove that for smaller values of the pressure-viscosity 
coefficient the traditional modification of the Reynolds equation is a very accurate approximation of the 
pressure field. 

It has been argued that the behavior of a piezoviscous fluid cannot be adequately described by a Reynolds 
type equation if the principal shear stress r is not less than the reciprocal of the pressure-viscosity coefficient 
a, cf. [18]. This is not surprising since problems of nonexistence and nonuniqueness are expected for the 
full balance of linear momentum equations if r > (a) -1 , cf. [7[. Similar conclusions can also be drawn from 
our modified Reynolds equations which ceases to be elliptic if the shear stress times the pressure-viscosity 
coefficient is larger than, or of the order of, one. 

In the piezoviscous regime, the pressure-dependence of viscosity can be modeled through the Barus 
relation (Barus [19]) 

H = Vo e ap , (1) 

where Ho denotes the dynamic viscosity measured at the ambient pressure (p = 0) and a is a positive 
parameter related to the rate of change of viscosity with respect to pressure, often referred to as the pressure- 
viscosity coefficient. The Barus relationship, together with the Roelands formula (Roelands [20]) 


H = Mo 



i—(i + v / pr) z 


( 2 ) 


where hr = 6.31 -10 5 Pa s, pu = 1.98 • 10 8 Pa and Z is a dimensionless parameter usually adjusted at 
ambient pressure to the Barus relation through the equation 


Z = 


a PR 

In pa - In HR ’ 


( 3 ) 


are the most widely used models in elastohydrodynamic lubrication. For applications of these and other 
experimentally validated formulas for the variation of viscosity with pressure, temperature and density, see, 
e.g., [2, 21, 22, 23, 3, 4], 

In general, the pressure-viscosity coefficient a depends on the lubricant, and on the pressure, temperature 
and shear rate in the contact area, cf. [3] . For different lubricant oils, its value has been shown to vary between 
1 • 10 -8 (Pa) -1 and 4-10 -8 (Pa) -1 , cf. [24, 25, 26]. If the constant Z in the Roelands formula (2) is computed 
from equation (3), then both the Barus and the Roelands relation give = a. We stress that although 
we rely on the assumption that the viscosity h varies with pressure p according to the Barus or Roelands 
formula, our computations hold for more general pressure-viscosity relationships for which is of the 
order of e or smaller. 

Before concluding this introductory section, it is worthwhile recognizing a marked departure of the 
constitutive relation of a piezoviscous fluid from that of the classical incompressible Navier-Stokes fluid 
with constant viscosity. While the latter is described by an explicit constitutive function for the stress in 
terms of the symmetric part of the velocity gradient, the former is an implicit relationship for the Cauchy 
stress and the symmetric part of the velocity gradient, leading to a totally different structure to the equations 
governing the flows of such fluids which in turn raises interesting issues with regard to both the mathematical 
and numerical analysis of the governing equations. 

A one-dimensional rate-type implicit model to describe the non-Newtonian response of fluids was intro¬ 
duced by Burgers [27]. His model includes as a special case the pioneering model to describe the viscoelastic 
response of fluids that was advanced by Maxwell [28]. Maxwell’s model is however not an implicit model 
as the symmetric part of the velocity gradient can be expressed explicitly in terms of the stress and the 
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time rate of the stress. Oldroyd [29] developed a systematic procedure to generate properly invariant three- 
dimensional rate type implicit constitutive relations. Such fluid models can be used to describe the flow of 
viscoelastic fluids and it is our aim to generalize the type of approximation that is being carried out here to 
include rate type implicit constitutive relations. 


2. Lubrication approximation 

Consider the following equations governing the isothermal flow of an incompressible, homogeneous, vis¬ 
cous fluid 

P + (v'V)vj-2V'( p{p) D(v)) + Vp = ph , (4) 


V • v = 0, (5) 

where p > 0 is the constant density of the fluid, v = {u, v, w) is the velocity field and p is a scalar variable, 
often referred to as the mechanical pressure, associated with the incompressibility constraint (5). Moreover, 
we have assumed above that the stress response is linear in the symmetric part of the velocity gradient, 
D(v) = \ (Vv + (Vv) T ), with the viscosity p depending (in isothermal conditions) only on the pressure p, 
so that the stress tensor T reads as 

T = —pi + 2p(p) D(v) . (6) 

Ignoring the body force b and limiting oneself to steady motions, equations (4)-(5) become 

p (v • V) v — p{p) Av — 2D(v) Vp(p) + Vp = 0, (7) 

V • v = 0 . (8) 

Note that (6) defines an implicit relationship f(T,D) = 0 between the stress tensor T and the symmetric 
part of the velocity gradient D since in an incompressible fluid trD = 0 so that p = — |trT, i.e. the 
mechanical pressure is just the mean normal stress. 

Let us assume, for expediency, that the viscosity depends on the pressure through the Barus relation 
(1). Introducing the non-dimensionalized (starred) quantities 

x* = L _1 x, v* = t/~ 1 v, p* = py 1 p, p* = P~ 1 p 1 a* = P a , 

where L and U represent typical length and velocity scales and the characteristic pressure is taken to be 
P = poUL -1 , we can rewrite equations (7)-(8) as 

p* (p*)A*v* + 2 a*p*ip*) D*(v*) V*p* - V*p* = Re (v* • V*)v*, (9) 


V* -v* =0, 


( 10 ) 


where Re = pULp g 1 denotes the usual Reynolds number. Note that, had we defined a more generally 


through 


1 dp 
a = - — 
p dp 


( 11 ) 


the form of the nondimensional system (9)-(10) would still remain the same. 
Let us restrict our attention to two-dimensional plane flows 


*(x*) = iu*ix*,y*),v*ix*,y*)) , p*=p*ix*,y*).. 
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so that equations (9)-(10) can be recast as 


'd 2 u* d 2 u* \ * * f du* dp* / du* dv*\ dp*\ dp* 

^ \dx* 2 dy* 2 J a ^ dx* dx* \dy* + dx* J dy* J dx* 


„ , * du* * du* 

= Re I u + d ^ , , 
dx* dy*J 


d 2 v* d 2 i 


dx 


*2 


dy 


*2 


* * 
a p 


du* 

dy* 


dv* 

dx* 


dp* 

dx* 


dv* dp* 
dy* dy* 


dp* 

dy* 


= Re ( u 


,dv* 

dx* 


,dv* 

dy* 


( 12 ) 


(13) 


du* dv* 
dx* + dy* 

We make the following assumptions 


(14) 


• the flow takes place between two almost parallel surfaces situated at y* = 0 and y* = H*(x*)\ 

• curvature of the surfaces can be neglected; 

• the lubrication film is thin, that is, H*{x*) = eh*(x*), where e <C 1 denotes a small non-dimensional 
parameter; 

• the characteristic lengths in x- and y-directions scale as L y /L x = 0(e); 

• the flow is slow enough or the viscosity high enough so that Re = 0(e); 

• the scaled viscosity is of the order of one, that is yi = 0(1). 


Remark 1. Choosing Re = 0(e) we end up neglecting the inertial effects from the outset, as usual in 
deriving the Reynolds approximation. At the same time, we simplify our presentation. In fact, given 
that the main pressure approximation is not affected by the inertial effects for Reynolds numbers up 
to the order 0(e -1 ), cf. [17] for computations in the constant-viscosity case, this simplification seems 
entirely warranted. 

As for the assumption n = 0(1), it is essential for the dimension reduction and, therefore, for the 
existence of a Reynolds type equation. It may not hold at the highest pressure regimes but, then again, 
for pressures higher than 0.5 GPa the entire viscosity-pressure relationships given by the Barus or 
Roelands formula becomes questionable, cf. [3, 25]. 
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Redefining the y-variable as y = e l y* and dropping the stars in (12)-(14), yields the system 


/ d 2 u _ 2 < 9 2 u\ f^dudp _ 2 c>u dp _ 1 dv dp\ 

\ dx 2 6 dy 2 J \ dx dx 6 dy dy + € dx dy) 


dp 

dx 


= Re 



-i 



2 ^u\ / dp du dp 2 dvdp\ _ -l^P 

\ dx 2 € dy 2 ) a ^ \ dy dx dx dx dy dy ) 6 dy 


= Re 




du 

■, dv 

— 

+ e“ — = 

dx 

dy 

c 

II 

d 

o 

, U° • n = 

v = U\ 

\] h n = 


0 , 

0, at y = 0, 

0, at y = h(x ), 


where U° and Uk denote the velocities of the parallel plates. We will write 


Re = e Re 


( 15 ) 


(16) 


(17) 

(18) 
(19) 


and assume that 


v(e,x) = e°v 1 (x) + e 1 v 2 (x) + e 2 v 3 (x) + ... , (20) 

p(e,x) = rV(x) + £ -V(x) + £ V(x) + ... , (21) 

where the functions v- 7 and p 7 and the parameter Re are of 0(1). 

There is still one dimensionless parameter, a , whose order has not been discussed. In the thinnest-film 
(elastohydrodynamic) lubrication problems, e is of the order ICR 6 or smaller, and the size of a should be 
taken to be of the order of e but in thicker film problems the order of the pressure-viscosity coefficient is e 2 
or smaller. We will consider these two cases separately. 

2.1. Case a = O(e) 

Substituting the asymptotic expressions (20)-(21) in equations (15)—(19) , writing 

a = eA , 
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where A = 0( 1), and keeping only the terms of the highest order, we obtain, at order e 2 in (15) and (16), 
at order e -1 in (17) and at order e° in (18) and (19) 


d 2 v} du 1 dp 1 dp 0 

+ A ^~^~ - = 0 ’ 

dy z dy dy dx 


( 22 ) 


d 2 v 1 /a^V 0 a ^ 1 ap 1 \ dp 1 n 

^ dy 2 + ^ \ dy dx + dy dy J dy 


(23) 


dv 1 

dy 


= 0 , 


(24) 


u 1 = U° , w 1 = 0, at 3 / = 0, 


(25) 


u 1 =U h , v 1 = 0, at y = h(x). 


Equations (24), (25) and (26), show that 

v x (x) =0. 

It then follows from (22) and (23) that 

d 2 ^ a 2 2 dp° (dv}\ 2 dp 0 

^ dy 2 ^ dx \ dy ) dx 


(26) 


(27) 


The next order equations read as 

d 2 u 2 dp 1 du 2 dp 1 du 1 dp 0 du 1 dp 2 

+ A ^~E~ ~E~ = ~E - 2A ^~^~ “b - 

dy z dy dy dx dx dx dy dy 


(28) 


d 2 v 2 „ , dp 1 dv 2 dp 2 „ a^ap 1 „ ap°a u 2 

+ 2 A ^^r A - ~a~ = ~ Afl ^r A - A ^~a~ 

dy z dy dy dy dy dx dx dy 


(29) 


dv 2 _ du 1 
dy dx ’ 


(30) 


u 2 = v 2 = 0 , at y = 0 , 


(31) 


v 2 = h'U h , u 2 = 0, at y = h(x ). (32) 

The one-dimensional Oseen problem (29)-(32) for (v 2 ,p 2 ) is solvable if and only if the compatibility condition 

r h du 1 


/ rrh 


dx 


d y = h U 


is satisfied. This condition can be written as 


r h 

j u 1 dy = 0 . 


d 

dx 


o 


6 


( 33 ) 







Multiplying (27) by y(y — h) and integrating across the film, yields 


r h 

2 n / u 1 dy + ph(U° 

Jo 


U h ) + A 2 y 2 


dp° 

dx 


y(y- h) 



h 3 dp 0 
6 dx 


Using (33) in (34) leads to the modified Reynolds equation 


d 

dx 



A 2 y / y(h-y) 



dp° 

dx 


"(U h -U°). 


2.2. Case a = 0(e 2 ) 

At the first order, we simply have 

dV _ dp° 
^ dy 2 dx 


( 34 ) 


(35) 


d 2 v 1 dp 1 

^ dy 2 dy ’ 


(36) 


dv 1 
dy 


= 0 , 


(37) 


= U°, 

v 1 = 0, 

at y = 0, 

(38) 

= u h , 

v 1 = 0, 

at y = h(x). 

(39) 


From (37)-(39) it again follows that tr(x) = 0. Equation (36) then yields dp 1 /dy = 0, i.e. the cross-film 
pressure also vanishes at this order. Therefore 


u 1 (x,y) = U° + 


U° — U h y(y — h ) dp° 
h ^ 2/i dx 


At the next order, we find that 


(40) 


d 2 u 2 dp 1 

^ dy 2 dx 


(41) 


d 2 v 2 dp 2 du 1 dp 0 

^ dy 2 dy ^ dy dx 


(42) 


dv 2 

du 1 


(43) 

dy 

dx 

5 

= v 2 

= 0 , 

at y = 0 , 

(44) 

u 2 

= 0 , 

at y = h{x). 

(45) 
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The one-dimensional Stokes problem (42)-(45) for (v 2 ,p 2 ) is solvable if and only if 


Substituting (40) into (46), yields 



du 1 
dx 


d y = h'U h . 


d_ 

dx 



(U° ~ U h ) ~2 


1 

2/a 


y{y - h ) d y 


dp° \ 

) 


+ tiu h = ti u h , 


which can be written as the classical Reynolds equation 


1 d / h 3 dp 0 \ 
12 da; \ p dx J 


j (U h + U°) . 


( 46 ) 


(47) 


At the next order, we obtain the following system for (v 3 ,p 3 ) 

d 2 u 3 dp 2 du 1 dp° du 1 dp 2 

^ dy 2 dx ^ dx dx ^ dy dy 


(48) 


d 2 v 3 dp 3 du 1 dp 1 

^ dy 2 dy ^ dy dx ’ 

dv 3 du 2 

dy dx ’ 

u 3 = v 3 = 0 , at y = 0 , 

u 3 = —^h' 2 U h , u 3 = 0 , at y = h(x). 

The one-dimensional Stokes problem (49)-(52) is solvable if and only if 


(49) 


(50) 

(51) 

(52) 



du 2 

dx 


dy = 0. 


On the other hand, from (41), (44) and (45) one concludes that 


u 2 (x,y) 


y(y ~ h{x )) dp 1 

2 /a dx 


This yields 

12 dx \ /i dx J 

We thus conclude that the Reynolds equation (47) provides a very accurate approximation for the pressure 
when a = 0(e 2 ) or smaller. 





3. Numerical results 


Contrary to [ 6 ], we do not simplify our modified Reynolds equation for numerical computations. Rather, 
we couple it with a nonlinear ODE governing the behavior of the along channel velocity and solve simulta¬ 
neously for the pressure and the main (along-channel) velocity component. The modified Reynolds equation 


d 

dec 



A 2 fi / y(h-y) 



j(U h -U°) 


(53) 


is a nonlinear second-order ODE for p°. Complemented with suitable boundary conditions, it becomes 
solvable when considered with equation 


d2yl , a 2 2 d P° 

M dy* + ' 4/i dx 1 

/du'\' 2 dp 0 

\dy ) ~ ~dx ’ 


(54) 


u 1 = U°, 

at y = 0 , 

(55) 


u 1 =U h , 

at y = h(x). 

(56) 


In order to illustrate the effect of the additional terms in the modified Reynolds equation, we study the 
classic model problem of a rigid cylinder rolling on a plane. (Refer to Szeri [4] for an overview of the problem 
statement and the traditional solution method.) Let h 0 be the minimum distance between the cylinder of 
radius R and the plane. Then the film thickness h(9) as a function of the angular coordinate 6 is given by 


h(6) = — — (1 + ncosd), n = —~———. 
w n y ’ h 0 + R 

The geometry of the problem is illustrated in Figure 1. 


(57) 



Figure 1: The geometry of the cylinder-plane assembly. 
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After performing the change of variables x = RsmO , the modified Reynolds equation (53) becomes 


d 

dO 


h 

12/z 


- a2 “L dv , 


dd d p° 


= nj h -U°) 


dh 

W 


where 


Similarly, for equation (54) we get 


1 


dd 


dx R cos 6 ’ 


^ = RsmO. 
dO 


dy 2 


2 2 dO dp 0 / du x \ dO dp 0 

^ dx dO \ dy ) dx dx 


Equation (58) will be subject to the Swift-Stieber boundary conditions (see Szeri [4]) 


0 


when 0 = — r, 


= 0 when 0 = 0 2) 


( 58 ) 


(59) 


(60) 


(61) 


where yet another unknown 0 2 is introduced through the unknown location of the exit boundary. To solve 
the resulting system of nonlinear equations, we perform iteratively the following three-step process 

Step 1. Compute (p°, 0 2 ) from (58), (61) with u 1 taken from the previous iteration (choose initially u 1 = 0). 
Step 2. Substitute p° into (60) and solve for u 1 on a set of predefined angular values within the interval 
[-j^TT, 021- 

Step 3. Interpolate u 1 between the chosen angular values and go to Step 1. 



Figure 2: Convergence of the difference of the pressure iterates in the L 2 -norm with a = 5.59 • 10 —8 (Pa) —1 and ho/R = 10 —4 . 
Here is the /c’th iteration in the procedure described in Steps 1—3 with p® corresponding to the solution of the modified 
Reynolds equation (58) when u 1 = 0. Given that the trailing edge angle 62 might change from one iteration to another, the 
L 2 -norm is computed using the current iterate’s O 2 and the previous pressure is taken as zero wherever the new 62 is larger 
than the previous one. 

The nonlinear equation (58) at Step 1 is approximated using a fixed-point iteration and the solution of 
the resulting linear equations is done by the finite element method with piecewise-linear elements. The finite 
element mesh is uniform in 0 and 10 4 elements were used. The correct value of 62 at which the condition 
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Figure 3: Convergence of the difference of the pressure iterates in the L 2 - norm when a = 1.75 • 10 ' (Pa) 1 and ho/R = 10 3 . 

(61) is approximately satisfied is sought through a binary search with initial lower and upper bounds at zero 
and respectively. 

The values of 0 at which problem (60), (55), (56) is solved in Step 2 are chosen to coincide with the nodes 
of the finite element mesh of Step 1. The nonlinear equation (60) is again approximated using a fixed-point 
scheme. The values of dp 0 /dd present in (60) are computed from the solution obtained at the previous 
step. The resulting linear equations are solved by the finite element method with 50 uniformly distributed 
piecewise-linear elements. Note that step (2) is perfectly parallelizable for different values of 9. 



6 [rad] 


Figure 4: The resulting pressure fields for the three models: standard Reynolds equation with a constant viscosity, standard 
Reynolds equation with a pressure-dependent viscosity and modified Reynolds equation with a pressure-dependent viscosity, 
the last two plotted with a = 5.59 • 10 -8 (Pa) -1 and ho/R = 10 -4 . Here O 2 ~ 0.0067rad and the maximum pressure values 
are 54.1 MPa and 53.5 MPa for the modified and standard Reynolds equations. 
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Figure 5: The resulting pressure fields for the three models: standard Reynolds equation with a constant viscosity, standard 
Reynolds equation with a pressure-dependent viscosity and modified Reynolds equation with a pressure-dependent viscosity, 
the last two plotted with a. = 1.75 • 10 —7 (Pa) -1 and ho/R = 10 —3 . Here 62 ~ 0.021 rad and the maximum pressure values are 
18.2 MPa for the modified and 15.9 MPa for the standard Reynolds equation. 



6 [rad] 


Figure 6: The viscosity distributions computed from the models with a pressure-dependent viscosity when a = 5.59-10 —8 (Pa) -1 
and ho/R = 10 —4 . The maximum viscosity values are 3.24Pas for the modified and 3.14Pas for the standard Reynolds 
equation. The corresponding pressure field is illustrated in Figure 4. 

A linear interpolation is performed in Step 3 after which the newly solved u 1 is substituted in (58). 
The convergence of the pressure between subsequent iterations in T-norm with the parameter values a = 
5.59 • 10 -8 (Pa)^ 1 and a = 1.75 • 10 -7 (Pa) -1 , ho/R = 10 -4 and ho/R = 10 3 , no = 0.158Pas, U° = Om/s 
and U h = lm/s are visualized in Figures 2 and 3. 

The resulting pressure and viscosity fields after nine iterations are compared to the solution of the 
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standard Reynolds equation with a pressure-dependent viscosity in Figures 4, 5, 6 and 7. For comparison, 
the isoviscous case is also depicted in Figures 4, 5. The results reveal that higher maximum pressure and 
viscosity values are obtained for both sets of parameter values when the modified Reynolds equation is 
applied. 



9 [rad] 


Figure 7: The viscosity distributions computed from the models with a pressure-dependent viscosity with a = 1.75* 10 -7 (Pa) -1 
and ho/R = 10 —3 . The maximum viscosity values are 3.8 Pas for the modified and 2.5 Pas for the standard Reynolds equations. 
The corresponding pressure field is illustrated in Figure 5. 

The pressure and viscosity differences predicted by the two models and shown in Figures 5 and 7 are 
larger than the ones in Figures 4 and 6. In Figures 5 and 7, a higher ratio /i 0 /R = 1CU 3 allows us to use 
a higher value of the pressure-viscosity coefficient which leads to lower pressures but, with our modified 
model, to a significantly higher viscosity. On the other hand, larger values of a, with a fixed h^/R ratio, 
caused numerical instabilities. The corresponding velocity fields u 1 are visualized in Figure 8. 

The computations were performed using MATLAB R2014a on a HP ProLiant BL460c server blade with 
two 8-core Intel Xeon E5-2670 processors and 256 GB of RAM provided by Aalto University IT Services. 


4. Conclusions 

We have proposed a new Reynolds type lubrication approximation for fluids with pressure dependent 
viscosities. Starting with the full balance of linear momentum equations, we have derived a coupled set 
of dimensionally reduced equations governing the flow of piezoviscous fluids in hydrodynamical lubrication 
problems. 

As shown, both with a rigorous analysis and through numerical computations, the correction to the 
classical Reynolds equation can be significant for certain values of the pressure-viscosity coefficient. The 
largest deviation from the classical lubrication approximation seems to occur for higher values of the pressure- 
viscosity coefficient and in less thin domains, at least if one adopts the Barus formula for the pressure-viscosity 
relationship. One should keep in mind though that the lower is the value of the pressure-viscosity coefficient 
and the thinner is the domain more concentrated, higher and less easily computable become the pressure 
spikes. 

Although derived for the Barus formula, our modified Reynolds approximation applies for more general 
pressure-viscosity relationship, e.g., the Roelands formula. Its usefulness beyond the study case presented 
here is under investigation. 
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6 [rad] 



6 [rad] 


Figure 8: The velocity field u 1 (0,y) corresponding to the pressure field of the final iteration with parameter values a = 
5.59 • 10 —8 (Pa) -1 , ho/R = 10 -4 (upper) and a = 1.75 • 10 -7 (Pa) -1 , ho/R = 10 —3 (lower). Note that the lower left area in 
both drawings corresponds to negative velocity suggesting backflow of lubricant before the pressure spike. 
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